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I INTRODUCTION 

In this lecture we will develop a concept of wakes and impedances for relativistic 
beams interacting with the surrounding environment. Among the numerous publi- 
cations and reviews on this subject, we refer here to recent books [1-3], where the 
reader can find a more detailed treatment and further references. 

We will use the CGS system of units throughout this paper. 

II INTERACTION OF MOVING CHARGES IN FREE 

SPACE 

We begin with interactions of particles that moving with constant velocity in 
free space. If the material walls are far from the particles, their effect in the first 
approximation can be neglected. 

Let us consider a leading particle of charge q moving with velocity v, and a 
trailing particle of unit charge moving behind the leading one on a parallel path at 
a distance s with an offset x, as shown in Fig. 1. We want to find the force which 
the leading particle exerts on the trailing one. 

We will use the following expressions for the electric and magnetic fields of a 
particle moving with a constant velocity (see, e.g., [4]): 

^=^§3. H=i.xi5, (1) 
where R is the vector drawn from point 1 to point 2, R* 2 = s 2 + x 2 /7 2 , and 

7 = (i-v 2 /c 2 y 1/2 . 

From Eq. (1) we find that the longitudinal force acting on the trailing charge is 

Fl = Ez = ~ 7 2( s 2 + 9 x 2/ 7 2)3/2 ' ( 2 ) 

and the transverse force is 
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FIGURE 1. A leading particle 1 and a trailing particle 2 traveling in free space with parallel 
velocities v. Shown also is the coordinate system x, z. 

In accelerator physics, the force F is often called the space-charge force. 

It is easy to see that for any position given by s and x, the longitudinal force 
decreases with the growth of 7 as 7 -2 . For the transverse force, if s 3> x/j, 
Ft ~ 7~ 4 , but for s = 0, F t ~ 7 . Hence, in the limit of ultrarelativistic particles 
moving parallel to each other, 7 — > 00, the electromagnetic interaction in free space 
vanishes. 

In this lecture, we will focus on the case of ultrarelativistic charges, where v — > c. 
The space-charge effects discussed above disappear in this limit, and the interaction 
between the particles is due only to the presence of material walls. 

Note that, taking the limit v — > c in Eq. (1) and recalling that s = vt — z, we 
can write the electromagnetic field of an ultrarelitivistic charge in free space as 

E = ^S(z -ct), H = z x E , (4) 

where r = xx + yy is a two-dimensional radius vector in a cylindrical coordinate 
system (x and y are the unit vectors in the directions of x and y, respectively). 

Ill PARTICLES MOVING IN A PERFECTLY 

CONDUCTING PIPE 

If particles from the above example move parallel to the axis in a perfectly 
conducting cylindrical pipe of arbitrary cross section, they induce image charges, 
on the surface of the wall, that screen the metal from the electromagnetic field 
of the particles. The image charges travel with the same velocity v (see Fig. 2). 
Since both the particles and the image charges move on parallel paths, in the limit 
v = c, according to the results in Section II, they do not interact with each other, 
no matter how close to the wall the particles are. 

Interaction between the particles in the ultrarelativistic limit can occur if 1) the 
wall is not perfectly conducting, or 2) the pipe is not cylindrical (which is usually 
due to the presence of RF cavities, flanges, bellows, beam position monitors, slots, 
etc., in the vacuum chamber). 




FIGURE 2. Particles traveling inside a perfectly conducting pipe of arbitrary cross section. 
Shown are the image charges on the wall generated by the leading charge. 

IV CAUSALITY AND THE "CATCH-UP" DISTANCE 



If a beam particle moves along a straight line with the speed of light, the elec- 
tromagnetic field of this particle scattered off the boundary discontinuities will not 
overtake it and, furthermore, will not affect the charges that travel ahead of it. The 
field can interact only with the trailing charges in the beam that move behind it. 
This constitutes the principle of causality in the theory of wakefields, according to 
which the interaction of a point charge moving with the speed of light propagates 
only downstream and never reaches the upstream part of the beam. 




FIGURE 3. A wall discontinuity located at z = scatters the electromagnetic field of an 
ultrarelativistic particle. When the particle moves to location z, the scattered field arrives to 
point z — s . 

We can estimate the distance at which the electromagnetic field produced by a 
leading charge reaches a trailing particles traveling at a distance s behind. Let us 
assume that a discontinuity located at the surface of a pipe of radius b at coordinate 
z = is passed by the leading particle at time t — 0, see Fig. 3. If the scattered 
field reaches point s at time t, then ct = yjjz — s) 2 + 6 2 , where z is a coordinate 
of the leading particle at time t, z = ct. Assuming that s < 6, from these two 



equations we find 



The distance z given by this equation is often called the catch-up distance. Only- 
after the leading charge has traveled this distance away from the discontinuity, can 
a particle at point s behind it feel the wakefield generated by the discontinuity. 



V ROUND PIPE WITH RESISTIVE WALLS 

Consider a round pipe of radius b, with finite wall conductivity a. A point charge 
moves along the z axis of the pipe with the speed of light, and a trailing particle 
follows the leading one at a distance s. Both particles are assumed to be on the axis 
of the pipe. Because of the symmetry of the problem, the only non-zero component 
of the electromagnetic field on the axis is E z , which, depending on the sign, either 
accelerates or decelerates the trailing charge. Our goal now is to find the field E z 
as a function of s. 

If the conductivity of the pipe is large enough, we can use perturbation theory 
to find the effect of the wall resistivity. In the first approximation, we consider the 
pipe as a perfectly conducting one. In this case, the electromagnetic field of the 
charge is the same as in free space and is given by Eqs. (4). For what follows, we 
will need only the magnetic field He, 

H e = —5{z-ct). (6) 

r 

Using the mathematical identity 

1 f°° 

5(z - ct) = — J dooe- iu){t - z/c \ (7) 
we will decompose Hq into a Fourier integral, 

/oo 
duH du) (r)e- mz/c , (8) 
-oo 



where 



HeM) = — ■ (9) 
nrc 



In the limit where the skin depth 5 corresponding to the frequency ou, 5 — 
cj y/2nauj, is much smaller than the pipe radius , 5 b, we can use the Leontovich 
boundary condition [5] that relates the tangential electric field E t on the wall with 
the magnetic one, 

E t = (Hx n, (10) 



where n is the unit vector normal to the surface and directed toward the metal, 
and 



Combining Eqs.(lO), (11) and (9), we find 



E zul \ r=b = -(l-i)J-^^-. (12) 

V o7T(T 7T0C 

Equation (12) gives us the longitudinal electric field on the wall, but we need the 
field on the axis of the pipe. To find the radial dependence of E zun we use Maxwell's 
equations, from which it follows that the electric field in a vacuum satisfies the wave 
equation. In the cylindrical coordinate system the wave equation for E z is 

1 d 2 E z (r,z,t) . , 

_ 1 8 2 E z {r, z, t) d 2 E z {r, z,t) Id ^(r, z, t) _ Q 
c 2 d 2 t d 2 z r dr Or 

Substituting the Fourier component E ZUJ (r)e~ tu; ^' z ^ c ^ into this equation, we find 

l_d_ r dE^_ = Q 

r dr dr 

This equation has a general solution E ZUI = A + B lnr, where A and are arbitrary 
constants. Since we do not expect E z to have a singularity on the axis, B = 0. 
Hence the electric field does not depend on r, E zu) = const, and 

E zu} \ r =0 = E zul \ r= b, (15) 

implying that E z J\ t= q is given by the same Eq. (12). Note that we have shown 
here that in the ultrarelativistic case the longitudinal electric field inside the pipe 
is constant throughout the pipe cross section. 

To find E z (z,t) we make the inverse Fourier transformation, 

/oo 
dwE^e-™ {t ~ z/c \ (16) 
-oo 

which gives 



E z (z,t) = (i - l)^r\U- I <L,yfce-*«->*\ (17) 
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The last integral can be taken analytically in the complex plane (see the Appendix), 
with the result 



for s > 0. For the points where s < 0, located in front of the charge, E z = in 
agreement with the causality principle. The positive sign of E z indicates that the 
trailing charge (if it has the same sign as q) will be accelerated in the wake. 

In our derivation we assumed that the magnetic field on the wall is the same as 
in the case of perfect conductivity. However, the magnetic field is generated not 
only by the beam current, but also by a displacement current, 

disp = J_^± (19) 

that vanishes in the limit of perfect conductivity. To be able to neglect the correc- 
tions to Hg due to j^ lsp , we must require the total displacement current to be much 
less then the beam current. In the Fourier representation, the time derivative d/dt 
reduces to multiplication by — icu, and the requirement is 

nb 2 ±-uE zu) = (20) 

or 



cj\3/2 Una 

7) «VlF' (21) 

In the space-time domain, the inverse wavenumber c/u corresponds to the distance 
s, and the condition of applicability of Eq. (18) is, 

/ cb 2 \ 1/3 

s » So =(— j . (22) 

The behavior of E z for very small values of s, s < sq, can be found in Ref. 6. Here 
we note only that the singularity in Eq. (18) saturates at small s, and the electric 
field changes sign and becomes negative at s = 0. This field decelerates the leading 
charge, as expected from the energy balance consideration. 



VI WAKE DEFINITION 

The electromagnetic interaction of charged particles in accelerators with the sur- 
rounding environment is usually a relatively small effect that can be considered 
as a perturbation. In the zeroth approximation, we can assume that the beam 
moves with a constant velocity along a straight line. We solve Maxwell's equation, 
find the fields, and then take into account the effect of these fields on a particle's 
motion. In this approach we neglect the second-order effects because the motion 
along the perturbed orbit can only slightly change the fields computed in the zeroth 
approximation. Those corrections are usually small, especially for ultrarelativistic 
particles. 



Another important feature of the interaction between the generated electromag- 
netic field and the particles is that in many cases of practical importance it is 
localized in a region small compared with the length of the beam orbit. It also 
occurs on a time scale much smaller than the characteristic oscillation times of the 
beam in the accelerator (such as betatron and synchrotron periods). This allows 
us to consider this interaction in the impulse approximation and characterize it by 
the amount of momentum transferred to the particle. 

Taking into account the above considerations, we will introduce the notion of the 
wake in the following way. Consider a leading particle 1 of charge q moving along 
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FIGURE 4. A leading particle 1 and a trailing particle 2 move parallel to each other in a vacuum 
chamber. 

axis z with a velocity close to the speed of light, v ~ c, so that z — ct (see Fig. 
4). A trailing particle 2 of unit charge moves parallel to the leading one, with the 
same velocity, at a distance s with offset p relative to the z-axis. The vector p 
is a two-dimensional vector perpendicular to the z-axis, p = (x,y). Although the 
two particles move in a vacuum, there are material boundaries in the problem that 
scatter the electromagnetic field and result in interaction between the particles. 

Let us assume that we solved Maxwell's equation and found the electromagnetic 
field generated by the first particle. We calculate the change of the momentum 
Ap of the second particle caused by this field as a function of the offset p and the 
distance s, 

/oo 
dt [E(p, z,t) + zx B(p, z, t)] z=ct _ s . (23) 
-oo 

Note that we integrate here along a straight line — the unperturbed orbit of the 
second particle. The integration limits are extended from minus to plus infinity, 
assuming that the integral converges. 

Since the beam dynamics is different in the longitudinal and transverse direc- 
tions, it is useful to separate the longitudinal momentum Ap z from the transverse 
component Ap±. With the proper sign and the normalization factor c/q, these two 



components are called the longitudinal and transverse wake functions (or simply 
wakes), 



c c f 

wi(p,s) = — Ap z = — / dtE z \ z=ct - s , 

w t (p, s) = ^Ap ± = ^Jdt[E ± + zx B] z=ct _ s . (24) 

Note the minus sign in the definition of wi — it is introduced so that the positive 
longitudinal wake corresponds to the energy loss of the trailing particle (if both 
the leading and trailing particles have the same sign of charge). The defined wakes 
have dimension cm -1 in CGS units and V/C in SI units. 1 

Because of the causality principle, the wakefield does not propagate in front of 
the leading charge, hence 

wi(p,s)=0, w t (p,s) = 0, fors<0. (25) 

It was assumed above that the electromagnetic field is localized in space and 
time and the integral in Eq. (23) converges. There are cases, however, where this 
is not true and the source of the wake is distributed uniformly along an extended 
path, such as the resistive wall wake of a long pipe, considered in Section V. In 
this case it is more convenient to introduce the wake per unit length of the path 
by dropping the integration in Eq. (23): 

w l{Pi s ) = --Ez\z=ct-s, 

q 

w t (p,s) = ^[E ± + zxB] z=ct _ s . (26) 

In this definition, the wakes acquire an additional dimension of inverse length, and 
has the dimension cm -2 in CGS and V/C/m in SI. 



VII PANOFSKY-WENZEL THEOREM 

Several general relations between longitudinal and transverse wakes can be ob- 
tained from Maxwell's equation without specifying the boundary condition for the 
fields. 

Let us introduce the vector R = (p, —s) (the negative sign in front of s is due 
to measuring s in the negative direction of z, see Fig. 4) and consider momentum 
Ap in Eq. (23) as a function of R. Let us assume that the electric and magnetic 
fields are specified through the vector potential A(r,t) and the scalar potential 

^ A useful relation between the units is 1 V/pC = 1.11 cm -1 . 



0(r, t), and compute Ap for the given fields. It is convenient to use the Lagrangian 
formulation of the equations of motion, 2 

ddL dL „ T 

= * s Vi - (27) 

with the Lagrangian for the trailing unit charge in the electromagnetic field 



v 2 1 

L = -mc 2 \jl - — + -Av-(j>. (28) 



Putting Eq. (28) into Eq. (27) yields (p = m7«) 

I (* + W = (29) 

Now, integrating this equation along the orbit of the trailing particle, x = const, 
y = const and z = ct — s, and assuming that the fields A and <fi vanish at infinity, 
we find 



Ap(R) = J 



dtV I -Av - 



^-V R W(R), (30) 

c 



where we introduced the wake potential W, 

W(R) = - fdt (-Av - 



q J \c 

= ^J dt{A z -<f>). (31) 

In the last equation we used v ~ cz. 

We just proved an important relation that states that all three components of 
the vector Ap can be obtained by differentiation of a single scalar function W . 
Recalling now the relation between the components of Ap and the wakes, Eq. 
(24), we find that 

dW dW 

= "ap? = IT- - = W < 32 > 



and hence 

dw t 



ds 



V P wi. (33) 



2 ^ This approach to the derivation of the Panofsky-Wcnzel theorem is due to A. Chao. 



This relation is usually referred to as the Panofsky-Wenzel theorem. Note that V p 
is a two-dimensional gradient with respect to coordinates x and y. 

One of the most important computational applications of the Panofsky-Wenzel 
theorem is that knowledge of the longitudinal wake function wi allows us to find 
the transverse wake w t by means of a simple integration of Eq. (33). 

We now prove another important property of W: it is a harmonic function of 
variables x and y, 

. „ T d 2 W d 2 W n 

s a? + -W = °- (34) 

To prove this, we will use the fact that both A and <fi satisfy the wave equation in 
free space, (d 2 /dt 2 - c 2 A)A = (d 2 /dt 2 - c 2 A)0 = 0. Hence 



c f , fd 2 



q J \dx 2 dy 2 J q J \dt 2 dz 2 

d 2 W d 2 W c f , ( d d\ ( d d 



/*(. 



+ - / ^ U: + Ut - c— - 0) . (35) 



dx 2 dy 2 q J \dt dz J \dt dz 
The last integral in this equation vanishes because 



and 



d d d „ d j . 

m +c Tz"m + vV= Jt (36) 



, /a a\ ( a a \ , , 

, d ( d d \ , . 
0. (37) 



It is interesting that the wake potential turns out to be a relativistic invariant. 
A covariant expression for it can be written as 

1 f°° 

W = — / A k u k dr, (38) 

Q J —CO 

where A k = (0, — A) is the 4- vector potential, u k = (cry, 071;) is the 4- vector velocity, 
and r is the proper time for the particle. 



VIII SYSTEMS WITH A SYMMETRY AXIS 



In Section VI we denned the wake as a function of the trailing particle offset 
relative to the path of the leading particle. In practical applications we are also 
interested in how the wake depends on the trajectory of the leading particle. We 
will assume that the system under consideration has a symmetry axis, and choose 
it as the z-axis of the coordinate system (see Fig. 5). Now the leading particle 1 
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FIGURE 5. Both the leading particle 1 and the trailing particle 2 are offset relative to the axis 
of the chamber. 

moves in the z direction with an offset given by vector p', and the trailing particle 
travels parallel to the leading one, with the same velocity, at a distance s behind 
the leading one, and with offset p relative to the axis. The vectors p' and p are 
the two-dimensional vectors perpendicular to the z-axis. The wake is still defined 
by Eq. (24), but now it will be considered as a function of p', p, and s 

wi = wi(p,p',s), 

w t = w t (p,p',s). (39) 

Usually the vacuum chamber is designed so that the system axis serves as an 
ideal orbit for the beam. Deviations from it are relatively small, and both vectors 
p and p' are typically much smaller than the size of the vacuum chamber. We can 
neglect them in wi and introduce the longitudinal wake function that depends only 
on s, 

wi(s)=wi(0,0,s). (40) 

If the vacuum chamber also has some symmetry elements (e.g., it has either 
circular, elliptical or rectangular cross section), the transverse wake on the axis, 
where p, p' = 0, vanishes, w t (0, 0, s) = 0. For small values of p, p' we can expand 
w t (p, p', s) keeping only the lowest-order linear terms. That gives a tensor relation 
between the transverse wake and the offsets, 



w t (p,p',s) = W 1 (s)p+W 2 (s)p', 



(41) 



where Wi and W2 are the two-dimensional tensors of rank 2. An example of 
the wake calculation for elliptical and rectangular cross sections of the pipe can be 
found in Ref. 7. 



IX AXISYMMETRIC SYSTEMS 

In an axisymmetric system W depends only on the absolute values of p, p', and 
the angle 9 between them. We can always chose a coordinate system such that the 
vector p' lies in the x-z plane (see Fig. 6), so that the potential function W will be 
a periodic even function of angle 9 in a cylindrical coordinate system. Decomposing 




FIGURE 6. Vectors p and p' in an axisymmetric system. 
W in Fourier series in 9 yields 



W(p,p',9,s) = J2 W m(p,p',s)cosm9. (42) 

m=0 

Putting this equation into Eq. (34) gives 

E(~^-^)cosme = 0, (43) 
^ \pdp dp p 2 J 

from which we can find an explicit dependence of W m of p, 

W m (p,p',s)=A m (p',s)p m . (44) 

In Eq. (44) we discarded a singular solution of Eq. (43) W m oc p~ m . 

It is also possible to find the dependence of W m versus p' (see [8]), which turns 
out to be 

A m (p',s)=F m (s)(p') m - (45) 
Using Eq. (32) we now find for the longitudinal and transverse wake functions 

wi = y^w} m) , w t = w\ m) (46) 



where 



wi m) = (pTp m FUs) cos m9, 

= m(p') m p m ~ 1 F m {.s) \rcosm6 - flsinmfll , (47) 

where r and are the unit vectors in the radial and azimuthal directions in the 
cylindrical coordinate system. Remember that in this equation we assume that the 
leading particle is in the plane 6 = 0. 

Equations (47) are valid for arbitrary values of p and p' . Near the axis, where the 
offsets are small, the higher-order terms with large values of m in these equations 
also become small. In this case, we can keep only the lower-order terms with m = 
(monopole) and m = 1 (dipole) wakes. For the monopole wake we find 

w l = w? S > = F' {8), (48) 

which shows that the longitudinal wake does not depend on the radius in an axisym- 
metric system. We also see that the monopole transverse wake vanishes, = 0. 

Since wf 1 ' does not depend on p', sometimes it is more convenient to compute 
the monopole wake for an offset orbit, p' ^ 0, rather than on the axis. 
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FIGURE 7. Vectors p and p 1 and unit vectors in cylindrical coordinate system. 

For the dipole wake (m = 1), the vector rcosd — 0sm9 lies in the direction of 
the x axis, that is in the direction of p', see Fig. 7. Hence, 

= p'F(s). (49) 

The wake given by Eq. (49) is usually normalized by the absolute value of the 
offset p', and the scalar function wf^/p' is called the transverse dipole wake w t , 

w t (s) = F{s). (50) 

Such a transverse wake has the dimension cm -2 or V/C/m. 3 In this definition, a 
positive transverse wake means a kick in the direction of the offset of the driving 
particle (if both particles have the same charge). 



3 ) If the original wake is denned per unit length, as in Eq. (26), then wt will have the dimension 
V/C/m 2 or cm -3 . 



X RESISTIVE WALL WAKE FUNCTIONS 



We are now in a position to find the wakes generated by a particle in a circular 
pipe with resistive walls. Using Eq. (18) and the definition Eq. (24) gives the 
longitudinal wake 



^-2^- <51) 

The minus sign here means that the trailing charge is accelerated in the wakefield. 

Let us now calculate the dipole transverse wake due to the resistive wall. First, 
we need to solve for the electromagnetic field of the leading charge q moving with 
an offset p' in a circular pipe. From the point of view of excitation of dipole modes, 
this charge can be considered as having a dipole moment d = qp'. In the zeroth 
approximation of perturbation theory, the electromagnetic field of a dipole moving 
with the speed of light in a perfectly conducting pipe is 



E = 25{z - ct) 



2{d-r)r-dr 2 d 
~4 + 7^ 



H = zxE. (52) 



The first term in the expression for E is a vacuum field of a relativistic dipole, and 
the second one is due to the image charges, which are generated in order to satisfy 
the boundary condition on the metal surface. 

Following the derivation in Section V, we find the magnetic field on the wall, 



Ad 

H e = — cos95(z-ct), (53) 
tr 



and take its Fourier transform, 



2d 

H eu , = —cos9, (54) 



where the angle 9 is measured from the direction of d. Then using the Leontovich 
boundary condition, Eq. (12), for the electric field E ZLU , 



^^ = -(1-0^^00861, (55) 
and making the inverse Fourier transformation, we obtain E z on the wall 



E z (z,p = bJ,t) = ^^ca*e, (56) 

where s = ct — z. Recalling that, according to Eq. (47) the dipole wake is a linear 
function of p, we conclude that 



E z (z,p, P ',t) = q -^J—co S 0^ (57) 



and the function F[(s) is 



which gives the following result for the transverse wake defined by Eq. (50): 

•*w = ;WI- (59) 

Analogous to the longitudinal wake, Eq. (18), this formula is valid only for s»so 
(see Eq. (22)). 



XI WAKEFIELD IN A BUNCH OF PARTICLES 

Up to now we have studied the interaction of two point charges traveling some 
distance s apart. If a beam consists of iV particles with the distribution function 
A(s) (defined so that X(s)ds gives the probability of finding a particle near point 
s), a given particle will interact with all other particles of the beam. To find the 
change of the longitudinal momentum of the particle at point s we need to sum the 
wakes from all other particles in the bunch, 

^Pz(s) = — ds'Xis'Ms' - s). (60) 

Here we use the causality principle and integrate only over the part of the bunch 
ahead of point s. In the ultrarelativistic limit the energy change AE(s) caused by 
the wakefield is equal to cAp z , so Eq. (60) can also be rewritten as 

/oo 
rfs , A(s')w / (s' -s). (61) 

Two integral characteristics of the strength of the wake are the average value of 
the energy loss AE av , and the rms spread in energy generated by the wake AE rms . 
These two quantities are defined by the following equations: 



AE av 



/oo 
dsAE(s)\(s), (62) 
-oo 



and 

1/2 



AE 



/oo 
ds(AE(s) - A£ av ) 2 A(s) 
-oo 



(63) 



As an example, let us calculate AE av and AE Tms for the resistive wall wake given 
by Eq. (51) and a Gaussian distribution function, 



s 2 



M = ^^\-**i> (64) 



where a s is the rms bunch length. Note that, since w in Eq. (51) is the wake per 
unit length of the pipe, we need to multiply the final answer by the pipe length L. 

A direction substitution of the wake Eq. (51) into Eq. (61) gives a divergent 
integral when s' s. 4 To remove the divergence, we need to recall that according 
to Eq. (48) the longitudinal wake is equal to the derivative of the longitudinal wake 
potential, wi = Fq(s) with F = (nb)^ 1 ^Jc/as for s > 0, and F = for s < 0. We 
then rewrite Eq. (61) as 



/oo 
ds'X(s') 
-oo 



.dF {s'-s) 



ds 



—Ne 2 L ds'^^-F (s' - s) 
Js ds 

Ne 2 Lc 1/2 ( s \ 

- -G - , (65) 
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where the function G(x) is 



G ^ = Wv^L ~V¥^T- (66) 

The plot of the function G(s/cr z ) is shown in Fig. 8, where the positive values 
of s correspond to the head of the bunch. We see that in the resistive wake the 
particles lose energy in the head of the bunch and get accelerated in the tail. On 
the average, of course, the losses overcome the gain. For the average energy loss 
one can find an analytical result: 

r(|) Ne 2 c 1/2 

^ = 23/%3/2 ba S/2 al/2 ■ ( 6? ) 

Numerical integration of Eq. (63) shows that the energy spread generated by the 
resistive wake is approximately equal to Ai£ av : 

A£ rms = 1.06A£ av . (68) 

If the beam is traveling in the pipe with an offset y relative to the axis, it will be 
deflected in the direction of the offset, by the transverse wakefields. To calculate 
the deflection angle 9 we use the relation 



4 ) The integral diverges because Eq. (51) is not valid for very small values of s, see Eq. (22). 
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FIGURE 8. Plot of the function G (s/a z ). 
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ds'\(s')w t (s f - s) 



NLr e yc 1/2 



-H - 



where the function H(x) is 

2 l/2 foo e -yV2 dy 



H(x) 



7T" 



3/2 Jx Vv- x 



The plot of the function H (s/a z ) is shown in Fig. 9. 

The deflection angle averaged over the distribution function is 



0„ 



£(i) NLr.yc 1 / 2 

21/2^/2^1/2^1/2' 



and the rms spread is 



((0- av )2)V2 = A l/2 



NLr e yc 1/2 



where 
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3\ r 2 (i/4)' 

4v^F _ 



(69) 
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(73) 



and if is the complete elliptic integral. The numerical value of A 1 / 2 is 0.186. 
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FIGURE 9. Plot of the function H (s/a z ). 



XII DEFINITION OF IMPEDANCE AND RELATION 



Knowledge of the longitudinal and transverse wake functions gives complete in- 
formation about the electromagnetic interaction of the beam with its environment. 
However, in many cases, especially in the study of beam instabilities, it is more con- 
venient to use the Fourier transform of the wake functions, or impedances. Also, it 
is often easier to calculate the impedance for a given geometry of the beam pipe, 
rather than the wake function. Recall, that in Section V we actually first com- 
puted the Fourier components of the wakes, and then, using the inverse Fourier 
transformation, found the wakes. 

For historical reasons the longitudinal Zi and transverse Z t impedances are de- 
fined as Fourier transforms of wakes with different factors, 



Note that the integration in Eqs. (74) can actually be extended into the region of 
negative values of s, because wi and Wt are equal to zero in that region. 

Impedance can also be defined for complex values of uo such that \muo > and 
the integrals, Eq. (74), converge. So defined, the impedance is an analytic function 
in the upper half-plane of the complex variable u>. 

We must keep in mind that other authors sometimes introduce definitions of the 
impedance that differ from the one given above. In Refs. 2 and 9 the longitudinal 



BETWEEN IMPEDANCE AND WAKE 




(74) 



impedance is denned as a complex conjugate to the one given by Eq. (74). Here 
we follow the definitions of Refs. 1 and 10. 

From the definitions in Eq. (74) it follows that the impedance satisfies the 
following symmetry conditions: 

ReZi(u) = ReZi(-u), lmZi(cu) = — lmZi(— u>), 

ReZ t {uj) = -ReZ t (-u), lmZ t (u) = lmZ t (-u). (75) 

The inverse Fourier transform relates the wakes to the impedances: 

i r°° 

w l (s) = — / duZ^e-™'* , 



Ms) = ^ I dwZ t {u)e-^ s ' c . (70) 



It turns out that the wakefield can actually be found if only the real part of the 
impedance is known. Indeed, we can rewrite Eq. (70) for wi as 

wi(s) — — / duj ReZi{uj) cos ImZ;(o;)sin— . (77) 

27T L c c J 



For negative values of s this formula should give wi = 0, hence 

/•OO 

= / du 

J oo 



' OO 

from which it follows that 

"OO 



„ / s WS u\s\ 
ReZi(u) cos h lmZi(u) sin ■ 



c c 



, , 1 f°° , n „ . , ws 2 f 00 , „ ^ . , us 
W[(s) = — / rfu;Rez;(ct;) cos — = — / duReZAu) cos — 
Joo c n J c 



(78) 



(79) 



A similar derivation for the transverse wake gives 

2 /"°° ujs 
Wt(s) — — / cL;ReZ t (u;) sin — . (80) 
k Jo c 

XIII ENERGY LOSS AND ReZ L 

We can relate the energy loss by the bunch to the real part of the longitudinal 
impedance. Indeed, 

/oo /»oo 
ds\(s) / ds'X{s') Wl (s' - s) 
-oo J —oo 

/OO poo i poo 

dsX(s) / ds'\(s')— / ^^(cu)e-^ (s '- s)/c 
-oo J —oo ^ " J —oo 

/oo 
c^(u;)|A(u;)| 2 , (81) 
-oo 
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where X(u) = ds\{s)e ius / c . Si 



ince A(— a;) = A* (a;), |A(cj)| 2 is an even function 




A£ av = / ^ReZ,(u;)|A(u;)| 2 , 



(82) 



where Q = Ne. 

For a point charge, A(s) = 5(s), A (a;) 



1, and the energy loss is 




L 



CO 



du;ReZ/(c<j). 



(83) 
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XIV KRAMERS-KRONIG RELATIONS 



Equations (79) and (80) relate KeZ(uj) to the wake function. Since Z(u) is given 
by Fourier transformation of w(s), the knowledge of ReZ(u) allows us to find Z(u), 
and hence ImZ(w). This means that haZ(u) and ReZ(u) are functionally related 
to each other. Mathematically this relation is manifested in the Kramers-Kronig 
dispersion relation, which can be written as 



where P.V. stands for the principal value of the integral. Taking the real and 
imaginary parts of this equation gives explicit relations between ReZ and ImZ: 



Assume that we have a solution of an electromagnetic problem corresponding 
to the current on the axis of a pipe with the time and space dependence given by 
e -iu)t+iujz/c^ Specifically, we know the electric field on the axis, E zu) (z)e~ tult . How 
can longitudinal impedance be found in terms of E ZUJ (z)l 

The longitudinal wake is equal to the integrated field E z generated by a point 
charge moving with the speed of light. The current corresponding to this point 
charge can be decomposed into a Fourier integral: 




(84) 




(85) 



XV USEFUL FORMULA FOR IMPEDANCE 

CALCULATION 




(86) 



Since we know the electric field generated by each harmonic, we can find the field 
due to the point charge as a superposition of E ZU} : 

E g {z,t) = ^-[ dujE ZLU (z)e-^. (87) 

^ J-oo 

For wi we then have 

i r°° 

M S ) = -7T- / dud(ct)E zw (ct - s)e- lwt 

^ J-oo 

= / dudzE zw (z)e-^ z+s)/c . (88) 

^ J-oo 

Comparing Eq. (88) with Eq. (74) we find 

POO 

Zt(u) = - / dzEUz)e- iuJz/c . (89) 



Hence the longitudinal impedance can be obtained simply by making Fourier trans- 
formation of E„Jz). 



XVI SMALL PILLBOX CAVITY IN ROUND PIPE 

As an example of using Eq. (89) we will derive here the longitudinal impedance 
Zi(uj) for a small axisymmetric cavity (pillbox) in a round perfectly conducting 
pipe, see Fig. 10. We assume that the wavelength associated with the frequency to 
is much larger than the dimension of the pillbox, c/u ^> g, h. 
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FIGURE 10. Small pillbox cavity in a round pipe. The dashed line near the wall shows the 
integration path in Eq. (92). 

First, we need to find the solution of Maxwell's equations corresponding to a unit 
current I u = e - tujt + tuJZ / c on the axis. Since the cavity is small, the magnitude of the 



magnetic field at the location of the cavity [z = 0) is approximately equal to H e 
on the wall of the pipe 



H, = 1. (90) 

Because E zu) does not depend on r (see Section V) we can choose the integration 
path in Eq. (89) close to the wall, as shown in Fig. 10, rather than the pipe axis. 
Along this path E ZUJ is not equal to zero only in the cavity gap, where \z\ ~ h, and 



/oo 
dzE zw (z)e 
-oo 



-ILOZ/C 



L 



dzEUz) = ~ = (91) 
c at c 



where <3> is the magnetic flux in the cross section of the cavity, $ = H e hg. As a 
result, 

Zfa) = -zu— = -*>Z —, (92) 

where Z = 4n/c = 377 Ohm. What we obtained is a purely inductive longitudinal 
impedance, which can be rewritten as 

Z,H = -i^C, (93) 
c z 

where the inductance C = 2gh/b. In CGS units the inductance has a dimension of 
cm, 1 cm = 1 nH. 

A more detailed calculation [11] shows that our method gives only an approximate 
solution of the problem. In addition to the solenoidal electric field generated by 
the time-varying magnetic flux in the cavity, there is also a contribution due to 
the potential component of the electric field. This results in a different numerical 
coefficient in Eq. (92) which depends on the ratio g/h. For example, for g = h the 
correction factor is 0.84. 



XVII INDUCTIVE IMPEDANCE 

We saw in the previous section that a small pillbox is characterized by inductive 
impedance if the frequency is not very large. This is a common feature of many 
small perturbations whose size is much smaller than the pipe radius (e.g., small 
holes, shallow obstacles on the wall, etc.) - for not very large frequencies their 
impedance is purely inductive. 



The longitudinal wake corresponding to the inductive impedance can be found 
by using Eq. (76) : 5 

w l (s)=C5'(s). (94) 

Because of the inductive wake, slices of the beam can change their energy, al- 
though the net energy loss for the bunch is zero because the real part of the 
impedance vanishes. We can find the energy change as a function of position 
within the bunch by using Eq. (61). Integration gives 

AE(s) = -e 2 NCX'{s). (95) 

For a Gaussian distribution function this reduces to 

e 2 NC je-?' 2 
°l V2ir 

where £ = s/a z . For the rms energy spread we find 



AE(s) = ~ ^= , (96) 



{AE 2 ) 1 ' 2 = 3" 3 / 4 (2tt)- 1 / 2 ^^. (97) 



o 2 



XVIII CAVITY IMPEDANCE 

In the more general case of a large cavity (Fig. 11) the beam excites cavity 
eigenmodes and the longitudinal wake in the cavity is composed of contributions 



FIGURE 11. An RF cavity with beam pipes. 

from single modes, 

to,(s) = £V n >(s). (98) 

n 



5 ' Since the integral involved in the calculation of wi(s) actually diverges at u — > oo, it should 
be treated as a generalized function. It is easier to verify Eq. (94) by putting it into Eq. (74) 
and checking that the resulting impedance is given by Eq. (93). 



For perfectly conducting walls, assuming that the modes do not propagate into the 
beam pipes, 6 the mode frequencies oo n are real. It should be no surprise that each 
partial wakefield oscillates with the frequency of the mode uo n , 



w 



(n) 



[s) = 2k n cos ^ 



00 n S 



(99) 



where k n is the loss factor, which depends on the geometry of the cavity and the 
mode number. As an example, for a cylindrical cavity with b = I and TMoio mode, 
^010 — 4.5/1. For a more rigorous derivation of the wake for a cavity, see Ref. 12. 

The cavity impedance for this wake can be calculated by using Eqs. (74). They 
give 



ReZ[ = nk n [S(u + oo n ) + 5{oo — oo n )], 
1 1 



ImZi = k n 



+ 



00 + 00 n 00 — 00 n 



(100) 



It is also easy to generalize the above wake for a cavity with lossy walls when 
the frequency of the mode has a small imaginary part j n , j n <C oo n . The wake now 
decays with time as 



wi(s) = 2A; n e- 7s/c cos ( 



00 n s 



;ioi) 



Again using Eq. (74), we can calculate the impedance. It has two peaks: one in 
the vicinity of oo = oo n and the other in the vicinity of oo = —oo n . Assuming that oo 
is close to oo n , we find 



Z, = 



7-i(u-w n y 



(102) 



REFERENCES 

1. A. W. Chao, Physics of Collective Beam Instabilities in High Energy Accelerators 
(Wiley, New York, 1993). 

2. B. W. Zotter and S. A. Kheifets, Impedances and Wakes in High-Energy Particle 
Accelerators (World Scientific, Singapore, 1998). 

3. A. W. Chao and M. Tigner, Handbook of Accelerator Physics and Engineering (World 
Scientific, Singapore, 1999). 

4. L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields, vol. 2 of Course of 
Theoretical Physics, 4th ed. (Pergamon, London, 1979) (Translated from the Rus- 
sian). 

5. L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, vol. 8 of 
Course of Theoretical Physics, 2nd ed. (Pergamon, London, 1960) (Translated from 
the Russian). 



6 ) This is true if the frequency of the mode is above the cutoff frequency for the pipe. 



6. K. L. F. Bane and M. Sands, The Short-Range Resistive Wall Wakefields, Tech. Rep. 
SLAC-PUB-95-7074, SLAC (December 1995). 

7. R. L. Gluckstern, J. van Zeijts, and B. Zotter, Phys. Rev. E47, 656 (1993). 

8. K. L. Bane, P. B. Wilson, and T. Weiland, in Proc. US Particle Accelerator School: 
Physics of Particle Accelerators, Upton, N.Y., 1983 (American Institute of Physics, 
New York, 1985), no. 127 in AIP Conference Proceedings, pp. 875-928. 

9. P. B. Wilson, in M. Month and M. Dienes, eds., Proc. US Particle Accelerator School: 
Physics of Particle Accelerators, Batavia, 1987 (American Institute of Physics, 
New York, 1989), no. 184 in AIP Conference Proceedings, pp. 525-564. 

10. S. A. Heifets and S. A. Kheifets, Review of Modern Physics 63, 631 (1991). 

11. S. S. Kurennoy and G. V. Stupakov, Particle Accelerators 45, 95 (1994). 

12. P. Wilson, High Energy Electron Linacs: Applications to Storage Ring RF Systems 
and Linear Colliders, Tech. Rep. SLAC-AP-2884 (Rev.), SLAC (November 1991). 



where £ — t — z/c — s/c. First, we change the integration variable uj to ( = uj^, 



We then consider ( as a complex variable. In order to make the integrand a single- 
valued function in the complex plane, we make a cut along the negative imaginary 
half-axis and deform the integration path from a straight line to a contour, shown 
in Fig. 12. The integral /(£) can now be split into two parts, h(C) and ^(O) 
corresponding to the left and right branches of the integration contour. 

In the first integral ii(£) we change the complex integration variable ( to the 
real positive variable r, ( — e 37U / 2 r = — it, so that — e^^y/r — (i — 1) \fr j \/2. 
This gives for l\ 



APPENDIX 



We show here how to calculate the integral in Eq. (17), 




(Al) 




(A2) 




2 3/2£3/2 



(A3) 



For the second integral ^(0 we choose the real positive variable r such that ( = 
e -m/2 r _ _^ r ^ w hi cn mea ns that ^/C — e~ m ^^Jr = (1 — i)y/r/y/2 and 




FIGURE 12. Complex plane £. The wavy line indicates the cut. The deformed integration 
contour consists of two paths, C\ and C2, corresponding to the integrals I\ and I2, respectively. 



= Hi). (A4) 

Hence 

m=uo+m) = -^^- (as) 



Substituting this equation into Eq. (17) gives Eq. (18). 



